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One  limit  to  the  fidelity  of  quantum  logic  operations  on  trapped  ions  arises  from  heating  of  the  ions’ 
collective  modes  of  motion.  Sympathetic  cooling  of  the  ions  during  the  logic  operations  may  eliminate  this 
source  of  errors.  We  discuss  the  benefits  and  drawbacks  of  this  proposal,  and  describe  possible  experimental 
implementations.  We  also  present  an  overview  of  trapped-ion  dynamics  in  this  scheme. 

PACS  number(s):  03.67. -a,  32.80.Pj 


I.  INTRODUCTION 

One  of  the  most  attractive  physical  systems  for  generating 
large  entangled  states  and  realizing  a  quantum  computer  [1] 
is  a  collection  of  cold  trapped  atomic  ions  [2].  The  ion  trap 
quantum  computer  stores  one  or  more  quantum  bits  (qubits) 
in  the  internal  states  of  each  trapped  ion,  and  quantum  logic 
gates  (implemented  by  interactions  with  externally  applied 
laser  beams)  can  couple  qubits  through  a  collective  quan¬ 
tized  mode  of  motion  of  the  ion  Coulomb  crystal.  Loss  of 
coherence  of  the  internal  states  of  trapped  ions  is  negligible 
under  proper  conditions  but  heating  of  the  motion  of  the  ion 
crystal  may  ultimately  limit  the  fidelity  of  logic  gates  of  this 
type.  In  fact,  such  heating  is  currently  a  limiting  factor  in  the 
National  Institute  of  Standards  and  Technology  (NIST)  ion- 
trap  quantum  logic  experiments  [3,4]. 

Electric  fields  from  the  environment  readily  couple  to  the 
motion  of  the  ions,  heating  the  ion  crystal  [3-7].  If  the  ion 
trap  is  much  larger  than  the  ion  crystal  size,  we  expect  these 
electric  fields  to  be  nearly  uniform  across  the  crystal.  Uni¬ 
form  fields  will  heat  only  modes  that  involve  center-of-mass 
motion  (COM  motion),  in  which  the  crystal  moves  as  a  rigid 
body.  Motional  modes  orthogonal  to  the  COM  motion,  for 
instance,  the  collective  breathing  mode,  require  field  gradi¬ 
ents  to  excite  their  motion.  The  heating  of  these  modes  is 
therefore  suppressed  [4].  However,  even  if  quantum  logic 
operations  use  such  a  “cold”  mode,  the  heating  of  the  COM 
motion  can  still  indirectly  limit  the  fidelity  of  logic  opera¬ 
tions.  Since  the  laser  coupling  of  an  internal  qubit  and  a 
motional  mode  depends  on  the  total  wave-packet  spread  of 
the  ion  containing  the  qubit,  the  thermal  COM  motion  can 
reduce  the  logic  fidelity  [3,4]. 

In  this  paper,  we  examine  sympathetic  cooling  [8]  in  a 
particular  scheme  for  which  we  can  continuously  laser  cool 
the  COM  motion  while  leaving  undisturbed  the  coherences 
of  both  the  internal  qubits  and  the  mode  used  for  quantum 
logic.  In  this  method,  one  applies  continuous  laser  cooling  to 
only  the  center  ion  of  a  Coulomb-coupled  string  of  an  odd 
number  of  ions.  One  can  address  the  center  ion  alone  if  the 
center  ion  is  of  a  different  ion  species  than  that  composing 
the  rest  of  the  string  [9].  Alternatively,  one  can  simply  focus 
the  cooling  beams  so  that  they  affect  only  the  center  ion.  In 
either  case,  the  cooling  affects  only  the  internal  states  of  the 


center  ion,  leaving  all  other  internal  coherences  intact.  If  the 
logic  operations  use  a  mode  in  which  the  center  ion  remains 
at  rest,  the  motional  coherences  in  that  mode  are  also  unaf¬ 
fected  by  the  cooling.  On  the  other  hand,  the  sympathetic 
cooling  keeps  the  COM  motion  cold,  reducing  the  thermal 
wave  packet  spread  of  the  ions.  In  the  following,  we  will 
discuss  the  dynamics  of  an  ion  string  in  which  all  ions  are 
identical  except  the  center  ion,  assuming  heating  by  a  uni¬ 
form  electric  field.  Our  results  give  guidelines  for  imple¬ 
menting  the  sympathetic  cooling  scheme.  Similar  results 
would  apply  to  two-  and  three-dimensional  ion  crystals  [10- 
13], 


II.  AXIAL  MODES  OF  MOTION 

We  consider  a  crystal  of  N  ions,  all  of  charge  q,  in  a  linear 
radiofrequency  (RF)  trap  [10,11].  The  linear  RF  trap  is  es¬ 
sentially  an  RF  quadrupole  mass  filter  with  a  static  confining 
potential  along  the  filter  axis  z.  If  the  radial  confinement  is 
sufficiently  strong  compared  to  the  axial  confinement,  the 
ions  will  line  up  along  the  z  axis  in  a  string  configuration 
[10,1 1].  There  is  no  RF  electric  field  along  z,  so  we  can  write 
the  axial  confining  potential  as  <p(z)  =  C\aoZ2/2  for  a0  a  con¬ 
stant.  The  potential  energy  of  the  string  is  then  given  by 


V(zi,  ■  ■  ■  ,zn ) 


,  N  2  N 
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for  Zi  the  position  of  the  ith  ion  in  the  string  (counting  from 
the  end  of  the  string).  The  first  term  in  the  potential  energy 
expresses  the  influence  of  the  static  confining  potential  along 
the  z  axis,  while  the  second  arises  from  the  mutual  Coulomb 
repulsion  of  the  ions.  For  a  single  ion  of  mass  m,  the  trap 
frequency  along  z  is  just  u>z=  sjc\a0/m. 

We  can  compute  the  equilibrium  positions  of  the  ions  in 
the  string  by  minimizing  the  potential  energy  of  Eq.  (1). 
Defining  a  length  scale  /  by  /'3  =  q/(47re0a0)  and  normal¬ 
izing  the  ion  positions  by  =  gives  a  set  of  equations 
for  the  as 
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which  has  analytic  solutions  only  up  to  N=3.  Steane  [1]  and 
James  [14]  have  computed  the  equilibrium  positions  of  ions 
in  strings  with  N  up  to  10.  The  potential  energy  is  indepen¬ 
dent  of  the  mass,  so  the  equilibrium  positions  of  ions  in  a 
string  are  independent  of  the  elemental  composition  of  the 
string  if  all  the  ions  have  the  same  charge. 

In  a  real  ion  trap  the  ions  will  have  some  nonzero  tem¬ 
perature  and  will  move  about  their  equilibrium  positions.  If 
the  ions  are  sufficiently  cold,  we  can  write  their  positions  as 
a  function  of  time  as  z,(f)=/ m,  +  qt(t),  where  qt(t)  is  small 
enough  to  allow  linearizing  all  forces.  We  focus  on  the  case 
of  an  odd  number  of  ions  N,  where  all  ions  have  mass  m, 
except  for  the  one  at  the  center  of  the  string  which  has  mass 
M.  The  ions  are  numbered  1  with  the  center  ion  la¬ 

beled  by  nc  =  (N+  l)/2.  Following  James  [14],  the  Lagrang- 
ian  for  the  resulting  small  oscillations  is 
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for  the  frequencies  £k  and  (orthonormal)  eigenvectors  v(k)  of 
the  N  normal  modes.  Because  of  our  normalization  of  the 
Lagrangian  (6),  the  £k  are  normalized  to  a>.  and  the  vik)  are 
expressed  in  terms  of  the  normalized  coordinates  2,(f).  In 
terms  of  the  physical  time  t,  the  frequency  of  the  Arth  mode  is 
£ka>z .  If  the  kth  mode  is  excited  with  an  amplitude  C,  we 
have 


qi(t)  =  Re[Cv\k)ei(S^,+  ^]  i*nc. 
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in  terms  of  the  physical  coordinates  qi(t). 

We  can  solve  for  the  normal  modes  analytically  for  N 
=  3 .  Exact  expressions  for  the  normal-mode  frequencies  are 
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normalized  to  wz .  The  mode  eigenvectors  are 
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We  define  a  normalized  time  as  T=coJ.  In  treating  the 
case  of  two  ion  species,  we  write  /x  =  M/m  for  the  mass  ratio 
of  the  two  species  and  normalize  the  amplitude  of  the  ion 
vibrations  q,(t)  as  Qi  =  qj\Jqa0,  ij=nc,  Qn  =q„c Vqflo/2-- 
The  Lagrangian  becomes 
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where 


AU  = 


Aij  i,j*nc 

Aij/'Jn  i  or  j  =  nc,i¥=j 


.Ajj/fi  i=j  =  nc 


(7) 


generalizing  the  result  of  James  [14]. 

The  Lagrangian  is  now  cast  in  the  canonical  form  for 
small  oscillations  in  the  coordinates  Qj(t).  To  find  the  nor¬ 
mal  modes,  we  solve  the  eigenvalue  equation 


u(1)  =  IV1(  l,^(13-5^),  l),  (14) 

u(2)=Al2(l,  0,  - 1),  (15) 

u(3)  =  JV3(l,^(13  — 5^),  l),  (16) 

in  terms  of  Qj(t).  Here  A, ,  /V2,  and  A3  are  normalization 
factors.  In  the  case  of  three  identical  ions  ( /x  =  1 ) ,  we  can 
express  the  mode  eigenvectors  in  terms  of  the  Q,(t)  as  v<r> 

=  (1,1,1)/V3,  u(2)  =  (l,0,-l)/\/2,  and  u(3)  =  (l,-2,l)/N/6. 
The  mode  eigenvectors,  in  this  special  case,  also  give  the  ion 
oscillation  amplitudes  in  terms  of  the  physical  coordinates 
q,(t).  Lor  three  identical  ions,  then,  pure  axial  COM  motion 
constitutes  a  normal  mode.  (This  result  holds  for  an  arbitrary 
number  of  identical  ions.)  We  also  note  that  the  center  ion 
does  not  move  in  mode  #2;  hence  the  frequency  and  eigen¬ 
vector  of  mode  #2  are  independent  of  /x.  Lor  any  odd  num¬ 
ber  N  of  ions  there  are  (N—  1  )/2  modes  for  which  the  center 
ion  does  not  move.  These  modes  will  likewise  have  frequen¬ 
cies  and  eigenvectors  independent  of  /x.  Moreover,  they  have 
v<n  >-m  =  ~v  + m  and  so  they  are  orthogonal  to  the  COM 
motion  and  do  not  couple  to  uniform  electric  fields.  The 
center  ion  moves  in  the  other  (A+l)/2  modes,  and  unless 
/x=l,  each  of  these  (A+l)/2  modes  has  a  component  of 
axial  COM  motion  and  therefore  couples  to  uniform  electric 
fields. 


032310-2 


SYMPATHETIC  COOLING  OF  TRAPPED  IONS  FOR  .  .  . 


PHYSICAL  REVIEW  A  61  032310 


k 


M- 


FIG.  1 .  Normalized  axial 
mode  frequencies  as  a  function  of 
fi  for  (a)  3,  (b)  5,  (c)  7,  and  (d)  9 
ions. 


For  N=  5  and  higher,  the  normal  mode  frequencies  de¬ 
pend  on  /r  in  a  complicated  way.  However,  it  is  easy  to  find 
the  frequencies  numerically.  Figure  1  shows  the  mode  fre¬ 
quencies  for  N=3,  5,  7,  and  9  as  a  function  of  //  for  0.01 
<//<  1 00.  The  modes  are  numbered  in  order  of  increasing 
frequency  (at  //  =  1 ) ,  and  are  normalized  to  coz .  In  each  case, 
the  lowest-lying  mode  has  all  ions  moving  in  the  same  di¬ 
rection  and  consists  of  pure  COM  motion  for  //  =  1 .  The 
even-numbered  modes  correspond  to  the  (N—  l)/2  modes 
for  which  the  center  ion  does  not  move.  Their  frequencies 
are  therefore  independent  of  fi.  For  both  very  large  and  very 
small  //  the  modes  pair  up,  as  shown  in  Fig.  1.  For  each  pair 
there  is  some  value  //>  1  for  which  the  modes  become  de¬ 
generate.  The  relative  spacing  between  modes  in  a  pair  is 
also  smaller  in  the  large-//  limit  than  in  the  small-//  limit. 

If  the  static  confining  potential  of  the  trap  is  not  perfectly 
harmonic,  the  normal  modes  of  motion  will  exchange  energy 
with  each  other.  This  problem  is  addressed  in  Ref.  [3]  (Sec. 
4.1.8).  These  effects  are  small  if  all  modes  are  reasonably 
cold.  Moreover,  mode  cross-coupling  is  a  resonant  process 
which  requires  the  mode  frequencies  to  be  related  as  sums  or 
differences.  Avoiding  such  resonances  by  tuning  the  trap  po¬ 
tentials  reduces  the  rate  of  energy  transfer  by  a  large  factor. 

The  Rabi  frequency  of  the  laser-ion  interaction  sets  the 
speed  of  quantum  logic  gates  and  the  linewidths  of  transi¬ 
tions  between  vibrational  states  of  the  ion  crystal.  If  we  per¬ 
form  quantum  logic  on  a  normal  mode  whose  frequency  is 
too  close  to  that  of  another  “spectator”  mode,  we  must  re¬ 
duce  our  Rabi  frequency,  and  thus  our  gate  speed,  to  avoid 
driving  transitions  on  the  spectator  mode;  otherwise  we  suf¬ 
fer  a  loss  of  fidelity.  Only  the  lowest-frequency  mode  is  well 
separated  from  all  other  modes  for  /t  very  large.  We  will 


show  below  (Sec.  IV)  that  this  mode  is  undesirable  for  use  in 
quantum  logic.  Hence,  in  order  to  maximize  gate  speed,  it  is 
best  to  use  a  cooling  ion  that  is  of  the  same  mass  as,  or 
lighter  than,  the  logic  ions.  In  this  case  mode  2  is  well  sepa¬ 
rated  from  all  other  modes,  as  shown  in  Fig.  1. 

III.  TRANSVERSE  MODES  OF  MOTION 

We  now  consider  the  motion  of  the  ions  transverse  to  the 
z  axis.  The  ions  experience  an  RF  potential  x  cos(( lt)(x2 
—y2)/ 2  for  a  suitable  choice  of  axes  x  and  y  perpendicular  to 
z,  where  17  is  the  frequency  of  the  RF  field  and  ^  is  a 
constant.  The  static  confining  potential  can  be  written 
(qa0/2)(z2—  ax2—  ( 1  —  a)y2)  at  the  position  of  the  ions 
(with  a  a  constant),  so  there  is  also  a  transverse  static  elec¬ 
tric  field.  To  analyze  the  ion  motion,  we  work  in  the  pseudo¬ 
potential  approximation  [15],  in  which  one  time  averages  the 
motion  over  a  period  of  the  RF  drive  to  find  the  ponderomo- 
tive  force  on  the  ion.  If  the  static  potential  is  negligible,  the 
RF  drive  gives  rise  to  an  effective  transverse  confining  po¬ 
tential  of  jmco2Q(x2  +  y2),  where  w,.0  =  q^/(  V2I7m)  for  an 
ion  of  mass  m.  If  we  include  the  effects  of  the  static  field,  the 
transverse  potential  becomes  km(u>2x2  + u>2y2),  where  cox 
l  —  aco~/cor0,  <j)y=  a>rQ yl  —  ( 1  —  a)a>z/a>~0.  Below 
we  will  assume  a=  1/2,  so  that  (oy=cox.  In  any  case,  the 
transverse  potential  is  that  of  a  simple  harmonic  oscillator,  as 
we  saw  also  for  the  axial  potential.  However,  the  transverse 
potential  depends  directly  on  the  ion’s  mass,  so  the  center 
ion  of  a  string  feels  a  different  trap  potential  than  the  others 
for  //  #  1 . 

We  define  e=  ,  so  that  cox=  u>zyje2—  1/2.  Then  the 
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normalized  Lagrangian  for  the  motion  along  x  is 


i  "  idxA2  i  J; 

‘iiy  (n) 

where  Xj=Xj\lc[a0  for  ii=nc  and  Xn  =xt  \jc\ct0/jL  are  normal¬ 
ized  ion  vibration  amplitudes  along  x.  Here 
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We  can  describe  the  normal  mode  frequencies  and  oscilla¬ 
tion  amplitudes  in  terms  of  the  eigenvectors  and  eigenvalues 
of  B'jj ,  just  as  for  the  axial  case  above.  The  normalizations  of 
the  time  and  position  coordinates  remain  the  same  as  in  the 
axial  case. 

In  the  previous  section,  we  assumed  that  the  radial  con¬ 
finement  of  the  ions  was  strong  enough  that  the  configuration 
of  ions  in  a  string  along  the  z  axis  was  always  stable.  How¬ 
ever,  for  sufficiently  small  e,  the  string  configuration  be¬ 
comes  unstable.  The  stable  configurations  for  different  val¬ 
ues  of  e  can  be  calculated  [16,17],  and  several  of  these 
configurations  have  been  observed  for  small  numbers  of  ions 
[10,11].  Rather  than  review  the  theory  of  these  configura¬ 
tions,  we  will  simply  find  the  range  of  validity  of  our  small- 
oscillation  Lagrangian  for  the  string  configuration.  The  string 
will  remain  stable  for  all  e  greater  than  some  es  =  es 

also  varies  with  N.  On  the  boundary  between  stable  and  un¬ 
stable  regions,  the  frequency  of  some  mode  goes  to  zero. 
Recalling  that  the  determinant  of  a  matrix  is  equal  to  the 
product  of  its  eigenvalues,  we  see  that  es(/r)  is  the  maxi¬ 
mum  value  of  e  satisfying  detB'(e,/r)  =  0  for  /r  fixed.  Fig¬ 
ure  2  shows  e5(/r)  as  a  function  of  /r  for  3,  5,  7,  and  9  ions. 
In  each  case,  there  is  a  cusp  in  e,(/r)  corresponding  to  the 
crossing  of  the  two  largest  solutions  to  detB'(e,yU,)  =  0.  The 
position  of  the  cusp  varies  with  the  number  of  ions,  but  lies 
between  /x  =  0. 1  and  /r  =  1  for  N^9.  The  positions  of  the 
cusps  are  labeled  with  arrows  in  Fig.  2.  For  /r  greater  than 
the  value  at  the  cusp,  e<es{/x)  corresponds  to  instability  of 
the  zigzag  mode,  so  that  the  string  breaks  into  a  configura¬ 
tion  in  which  each  ion  is  displaced  in  the  opposite  direction 
to  its  neighbors  [16,17].  For  /jl  smaller  than  the  value  at  the 


FIG.  2.  Trap  anisotropy  at  instability  of  the  string  configuration 
as  a  function  of  /r  for  3,  5,  7,  and  9  ions.  Arrows  indicate  the  cusps 
discussed  in  the  text. 

cusp,  es  is  independent  of  /jl.  In  this  regime,  e<  es  creates  an 
instability  in  a  mode  similar  to  the  zigzag  mode,  except  that 
the  center  ion  remains  fixed. 

We  can  proceed  to  calculate  the  frequencies  of  the  trans¬ 
verse  modes  for  values  e>  es(  j±).  Again,  these  frequencies 
are  normalized  to  the  axial  frequency  of  a  single  ion  of  mass 
m.  Figure  3  shows  the  transverse  mode  frequencies  for  3,  5, 
7,  and  9  ions  as  a  function  of  /r,  where  e  is  taken  equal  to 
l.les(/r).  The  modes  are  numbered  in  order  of  increasing 
frequency  at  fi  =  I  (all  ions  identical).  In  this  numbering 
scheme,  the  central  ion  moves  in  odd-numbered  modes  but 
not  in  even-numbered  modes.  The  frequencies  of  the  even- 
numbered  modes  appear  to  depend  on  /jl  because  they  are 
calculated  at  a  multiple  of  e5(/r);  for  constant  e  these  fre¬ 
quencies  are  independent  of  fi.  The  cusps  in  the  mode  fre¬ 
quencies  in  Fig.  3  arise  from  the  cusps  of  es( //)  at  the  cross¬ 
over  points  between  the  two  relevant  solutions  of  dct  B ' 
=  0.  Mode  frequencies  plotted  for  a  constant  value  of  e  do 
not  exhibit  these  cusps.  As  in  the  case  of  axial  motion,  the 
mode  frequencies  form  pairs  of  one  even-  and  one  odd- 
numbered  mode  for  small  /r.  However,  for  large  fi  all  but 
one  of  the  transverse  modes  become  degenerate.  The  only 
nondegenerate  transverse  mode  in  this  case  is  the  zigzag 
mode.  In  general,  the  modes  are  most  easily  resolved  from 
their  neighbors  for  /i=  1 ,  as  in  the  case  of  axial  motion. 
Increasing  e  reduces  the  frequency  spacing  between  nearly 
degenerate  modes.  At  e=l.leJ(/r)  and  /x=  1 ,  for  instance, 
the  fractional  spacing  between  the  cold  transverse  mode  of  3 
ions  and  its  nearest  neighbor  is  0.20,  but  for  e=  l.5es(fi)  the 
same  spacing  is  0.09. 

The  near  degeneracy  of  the  modes  for  large  or  small  /jl 
and  for  eles  significantly  greater  than  1  limits  the  usefulness 
of  these  modes  because  of  possible  mode  cross-coupling,  just 
as  for  the  axial  modes.  Resolving  a  particular  transverse 
mode  requires  operating  the  trap  near  the  point  at  which  the 
string  configuration  becomes  unstable,  i.e.,  e  near  e0(/j,).  In 
this  regime,  the  collective  motion  of  the  ions  is  quite  sensi- 
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FIG.  3.  Normalized  frequen¬ 
cies  of  the  transverse  modes  as  a 
function  of  /jl  with  e=l.le0(/r) 
for  (a)  3,  (b)  5,  (c)  7,  and  (d)  9 
ions. 


M- 


tive  to  uncontrolled  perturbations,  which  may  pose  signifi¬ 
cant  technical  problems  for  using  a  transverse  mode  in  quan¬ 
tum  logic  operations. 

IV.  MODE  HEATING 

Stochastic  electric  fields  present  on  the  ion  trap  elec¬ 
trodes,  for  instance,  from  fluctuating  surface  potentials,  can 
heat  the  various  normal  modes  of  motion  incoherently.  For 
ion  trap  characteristic  dimension  dtrap  much  larger  than  the 
size  of  the  ion  crystal  dions ,  these  fields  are  approximately 
uniform  across  the  ion  crystal,  so  they  couple  only  to  the 
COM  motion.  The  (N—  l)/2  even-numbered  modes  are  or¬ 
thogonal  to  the  COM  motion,  so  they  are  only  heated  by 
fluctuating  electric-field  gradients.  The  heating  rates  of  these 
modes  are  reduced  by  a  factor  of  at  least  (dionsldtrap)2<^  1 
as  compared  to  the  heating  of  the  other  modes  [4].  In  the 
following,  therefore,  we  will  neglect  the  effects  of  fluctuat¬ 
ing  field  gradients,  so  that  the  even-numbered  modes  do  not 
heat  at  all. 

The  analysis  of  Secs.  II  and  III  shows  that  the  motion  of  a 
crystal  of  N  ions  is  separable  into  the  3  N  normal  modes, 
each  of  which  is  equivalent  to  a  simple  harmonic  oscillator. 
Hence  we  can  quantize  the  crystal  motion  by  quantizing  the 
normal  modes.  The  Oh  normal  mode  gives  rise  to  a  ladder  of 
energy  levels  spaced  by  h£ku>z ,  with  3 N  such  ladders  in  all. 
If  we  now  write  the  uniform  electric-field  power  spectral 
density  as  SE(w),  we  can  generalize  the  result  of  Ref.  [18]  to 
give 

N 

+  2  vf 

j= i 

j*nc 


-  P  SE(Ck 

,l  k  4mh£ktoz 


/  vw 

nc 


for  the  heating  rate  of  the  Oh  mode,  expressed  in  terms  of 
the  average  number  of  quanta  gained  per  second.  Recall  that 
v\kl  is  the  oscillation  amplitude  of  the  z'th  ion  in  the  Ath 
normal  mode,  expressed  in  the  normalized  coordinates.  It  is 
useful  to  normalize  the  heating  rate  in  Eq.  (20)  to  the  heating 
rate  of  the  lowest-lying  axial  mode  of  a  string  of  identical 
ions.  This  normal  mode  consists  entirely  of  COM  motion 
and  we  write  v(jOM=  1  / \  A/  for  all  ions.  The  normalized  heat¬ 
ing  rate  of  the  Oh  mode  is  then 
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(21) 


where  we  have  assumed  that  the  spectral  density  SE(u>)  is 
constant  over  the  frequency  range  of  the  normal  modes,  i.e., 
SE(0Jz)  =  SE(£k(Dz). 

Figure  4  shows  plots  of  the  normalized  heating  rates  of 
the  axial  modes  for  N=3,  5,  7,  and  9  as  a  function  of  j±. 
Figure  5  is  the  same,  but  for  the  transverse  modes,  with  e 
=  1 . 1  ev .  The  numbering  of  modes  on  the  plots  of  heating 
rate  matches  the  numbering  on  the  corresponding  plots  of 
mode  frequency  (Figs  1  and  2). 

In  both  axial-mode  and  transverse-mode  plots,  the  even- 
numbered  modes  have  the  center  ion  at  rest,  while  the  center 
ion  moves  for  all  odd-numbered  modes.  We  see  from  Figs.  4 
and  5  that  the  modes  for  which  the  center  ion  is  fixed  can 
never  heat,  while  all  the  other  modes  always  heat  to  some 
extent  for  /z#  1 .  We  will  refer  to  these  modes  as  “cold”  and 
“hot”  modes,  respectively.  If  the  ions  are  identical,  only  the 
modes  with  all  ions  moving  with  the  same  amplitude  (COM 
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FIG.  4.  Normalized  heating 
rates  for  the  axial  modes  as  a 
function  of  /x  for  (a)  3,  (b)  5,  (c) 
7,  and  (d)  9  ions. 


modes)  can  heat.  There  are  three  such  modes,  one  along  x, 
one  along  y,  and  one  along  z.  In  interpreting  Figs.  4  and  5,  it 
is  important  to  recall  that  the  normalized  heating  rate  defined 
in  Eq.  (21)  is  inversely  proportional  to  the  mode  frequency. 
For  instance,  the  /jl  dependence  of  the  heating  rate  of  the 
highest-frequency  transverse  mode  can  be  largely  ascribed  to 
variations  in  the  mode  frequency,  rather  than  to  changes  in 
the  coupling  of  the  mode  to  the  electric  field. 

V.  PROSPECTS  FOR  SYMPATHETIC  COOLING 

Heating  reduces  logic  gate  fidelity  in  two  ways.  The  logic 
mode  itself  can  be  heated,  but  by  choosing  a  cold  mode,  we 
can  render  this  effect  negligible.  On  the  other  hand,  the  Rabi 
frequency  of  the  transition  between  logic-mode  motional 
states  depends  on  the  total  wave-packet  spread  of  the  ion 
involved  in  the  transition  [3,4].  Heating  on  modes  other  than 
the  logic  mode  can  thus  lead  to  unknown,  uncontrolled 
changes  in  this  Rabi  frequency,  resulting  in  overdriving  or 
underdriving  of  the  transition.  The  purpose  of  sympathetic 
cooling  is  to  remove  this  effect  by  cooling  the  center  ion  and 
thus  all  hot  modes. 

In  the  foregoing,  we  have  chosen  to  consider  only  the 
case  of  a  crystal  of  an  odd  number  of  ions,  with  the  cooling 
ion  at  the  center.  We  now  see  that  this  is  the  only  case 
suitable  for  sympathetic  cooling,  since  only  in  this  case  do 
we  find  both  (a)  cold  modes  for  arbitrary  /x  and  (b)  isolation 
of  motion  of  the  cold  modes  from  motion  of  the  cooling  ion. 
As  long  as  the  crystal  is  symmetric  under  inversion  in  z,  the 
normal  modes  must  be  either  symmetric  (hot)  or  antisym¬ 
metric  (cold)  under  inversion  in  z  .  If  the  cooling  ion  is  not  at 
the  center  of  the  crystal,  the  crystal  symmetry  is  broken  and 


cold  modes  will  only  appear  for  particular  values  of  /x.  The 
case  of  a  crystal  of  an  even  number  of  ions,  with  two  cooling 
ions  at  the  center,  will  again  yield  cold  modes.  However, 
excitation  of  any  mode  will  cause  the  cooling  ions  to  move, 
so  that  the  cold  modes  are  not  well  isolated  from  the  sympa¬ 
thetic  cooling  in  this  case. 

For  sympathetic  cooling  to  be  useful,  we  must  find  a  cold 
mode  suitable  for  use  in  quantum  logic.  The  cold  mode  must 
be  spectrally  well  separated  from  any  other  modes  in  order  to 
maximize  gate  speed.  We  can  use  the  lowest-lying  cold  axial 
mode  as  the  logic  mode  for  /z;£3.  In  this  mode,  called  the 
breathing  mode,  the  center  ion  remains  fixed  and  the  spac- 
ings  between  ions  expand  and  contract  in  unison.  Unless  the 
trap  is  operated  very  close  to  the  instability  point  of  the 
string  configuration,  the  breathing  mode  is  better  separated 
from  its  neighbors  than  are  any  of  the  cold  transverse  modes. 
For  /i£3  any  cold  mode,  either  axial  or  transverse,  is  nearly 
degenerate  with  a  hot  mode.  In  this  regime  one  must  make  a 
specific  calculation  of  mode  frequencies  in  order  to  find  the 
best-resolved  cold  mode.  Even  so,  the  cold  axial  modes  are 
again  better  separated  from  their  neighbors  than  are  the  cold 
transverse  modes,  except  for  e  very  close  to  es{  /x).  It  seems 
best  to  select  a  cold  axial  mode  as  the  logic  mode  in  most 
cases. 

By  selecting  our  laser-beam  geometry  appropriately,  we 
can  ensure  that  the  Rabi  frequency  of  the  motional  transition 
on  the  axial  mode  used  for  logic  depends  chiefly  on  the 
spread  of  the  ion  wave  packet  along  z.  In  this  case,  heating  of 
the  axial  modes  will  affect  logic-gate  fidelity,  but  heating  of 
the  transverse  modes  will  have  little  effect.  If  the  mass  of  the 
central  ion  is  nearly  the  same  as  that  of  the  others  ( /x  1 ) , 
only  the  lowest  axial  mode  will  heat  significantly,  and  we 
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FIG.  5.  Normalized  heating 
rates  for  the  transverse  modes  as  a 
function  of  /x  with  e=l.le0(/z) 
for  (a)  3,  (b)  5,  (c)  7,  and  (d)  9 
ions. 


can  continuously  cool  this  mode  by  cooling  only  the  central 
ion,  ensuring  that  all  ions  remain  in  the  Lamb-Dicke  limit 
[3].  If  /z  is  not  near  1,  we  must  cool  all  (N+  l)/2  hot  modes 
(again  by  addressing  the  central  ion)  to  keep  all  ions  in  the 
Lamb-Dicke  limit. 

The  analysis  above  indicates  that,  all  other  things  being 
equal,  we  are  best  off  if  our  substituted  ion  is  identical  to,  or 
is  an  isotope  of,  the  logic  ions.  However,  sympathetic  cool¬ 
ing  can  still  be  useful  if  the  two  ion  species  have  different 
masses.  For  example,  we  can  consider  sympathetic  cooling 
using  the  species  9Be+  and  24Mg+.  Linear  traps  constructed 
at  NIST  have  demonstrated  axial  secular  frequencies  of  over 
10  MHz  for  single  trapped  9Be+  ions.  For  three  ions  with 
24Mg+  as  the  central  ion,  oj_(Be  1 )  =  2  ttX  10  MHz  yields  a 
spacing  of  1.6  MHz  between  the  cold  axial  breathing  mode 
and  its  nearest  neighbor.  If  we  reverse  the  roles  of  the  ions 
[w;.(Mg+)  =  27rX  10  MHz],  the  spacing  increases  to  6.2 
MHz.  The  transverse  modes  are  much  harder  to  resolve  from 
each  other.  For  three  ions  with  24 3 4Mg+  in  the  center,  we 
require  a>,.0(Be+)  =  2 -77X27.6  MHz  to  obtain  e=  1 . 1  es ,  and 
the  spacing  between  the  cold  transverse  zigzag  mode  and  its 
nearest  neighbor  is  only  560  kHz.  Reversing  the  roles  of  the 
ions,  we  find  6=1.16^  at  wr0(Mg+)  =  2ttX  14.7  MHz  with 
a  spacing  of  1.1  MHz.  For  this  combination  of  ion  species, 
the  cold  axial  breathing  mode  seems  most  appropriate  for 


logic.  For  a  string  of  3  or  5  ions,  sympathetic  cooling  would 
require  driving  transitions  on  2  or  3  axial-mode  sidebands, 
respectively.  From  this  example  we  see  that  sympathetic 
cooling  can  be  useful  even  for  ion  mass  ratios  of  nearly  3 
to  1. 

VI.  CONCLUSION 

We  have  investigated  a  particular  sympathetic  cooling 
scheme  for  the  case  of  an  ion  string  confined  in  a  linear  RF 
trap.  We  have  numerically  calculated  the  mode  frequencies 
of  the  axial  and  transverse  modes  as  functions  of  the  mass 
ratio  /x  and  trap  anisotropy  e  for  3,  5,  7,  and  9  ions.  We  have 
also  calculated  the  heating  rates  of  these  modes  relative  to 
the  heating  rate  of  a  single  ion,  assuming  that  the  heating  is 
driven  by  a  uniform  stochastic  electric  field.  The  results  in¬ 
dicate  that  the  scheme  is  feasible  for  many  choices  of  ion 
species  if  we  use  a  cold  axial  mode  as  the  logic  mode.  The 
optimal  implementation  of  the  scheme  employs  two  ion  spe¬ 
cies  of  nearly  equal  mass.  However,  a  demonstration  of  sym¬ 
pathetic  cooling  using  9Be+  and  24Mg+  appears  well  within 
the  reach  of  current  experimental  technique. 
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